Hydrogen permitted lines in the first near-IR spectra of Th 28 
microjet: accretion or ejection tracers? 

Deirdre Coffey^, Francesca Bacciotti^, Linda Podio^, Brunella Nisini^ 



ABSTRACT 

We report the first near-infrared detection of tiie bipolar microjet from T Tauri 
star TtiA 15-28 (hereafter Th 28). Spectra were obtained with Very Large Tele- 
scope(VLT) /ISAAC for the slit both perpendicular and parallel to the flow to examine 
jet kinematics and gas physics within the flrst arcsecond from the star. The jet was 
successfully detected in both molecular and atomic lines. The H2 component was found 
to be entirely blueshifted around the base of the bipolar jet. It shows that only the blue 
lobe is emitting in H2 while light is scattered in the direction of the red lobe, highlight- 
ing an asymmetric extinction and/or excitation between the two lobes. Consistent with 
this view, the red lobe is brighter in all atomic lines. Interestingly, the jet was detected 
not only in [Fe II], but also in Br7 and Pa/3 lines. Though considered tracers mainly of 
accretion, we find that these high excitation hydrogen permitted lines trace the jet as 
far as 150 AU from the star. This is confirmed in a number of ways: the presence of 
the [Fe II] 2.13 /um line which is of similarly high excitation; H I velocities which match 
the jet [Fe II] velocities in both the blue and red lobe; and high electron density close 
to the source of >6xl0^ cm~^ derived from the [Fe II] 1.64,1.60 fim ratio. These near- 
infrared data complement Hubble Space Telescope Imaging Spectrograph (HST/STIS) 
optical and near- ultraviolet data for the same target which were used in a jet rotation 
study, although no rotation signature could be identified here due to insufficient angu- 
lar resolution. The unpublished HST/STIS Ha emission is included here along side the 
other H I lines. Identifying Br7 and Pa/3 as tracers of ejection is significant because of 
the importance of finding strong near-infrared probes close to the star, where forbidden 
lines are quenched, which will help understand accretion-ejection when observed with 
high spatial resolution instruments such as VLTI/ AMBER. 

Subject headings: ISM: jets and outflows — stars: formation, — stars: individual (ThA 15-28) 
— stars: pre-main sequence 



^The Dublin Institute for Advanced Studies, 
Dublin 2, Ireland 

^Osservatorio Astrofisico di Arcetri, Largo E. Fermi 
5, 50125 Firenze, Italy 

^Osservatorio di Roma, Frascati, Roma, Italy 



1. Introduction 

Both observational and theoretical studies 
of the star formation process have unveiled 
a complex picture of simultaneous accretion 
and ejection which is yet to be fully under- 
stood. It is clear, however, that there ex- 



ists a tight correlation between the two mech- 
anisn is, as predicted b y theoretical models 
(e.g. , iKonigl et a l."200d'l and support ed by ob- 
servations (e.g.. Harti gan et al.lll995l ). 

Explaining the underlying physics rests ul- 
timately on identifying the origin of observed 
emission of different atomic species, in or- 
der to determine the gas properties, structure 
and dynamics. For example, the diffuse gas in 
protostellar outflows produces numerous per- 
mitted and forbidden lines, whose excitation 
properties are reasonably well known. Obser- 
vations of forbidden emission lines (such as 
0°, N"*" and S"^ at optical wavelengths and 
Fe"*" in near-infrared (near-IR) wavelengths) 
show that they trace hig h velocity, low den - 
sity outflow activity (e.g., Podio et al. 20061 ) . 
Permitted hydrogen lines, however, are more 
complex to analyse and so their exact ori- 
gin is unclear. Advancements in understand- 
ing often rely on correlations and modeling 
rather than direct observation, particularly 
concerning the sub-arcsecond spatial scale on 
which the accretion/ejection engine operates. 
Hydrogen line profiles have been observed 
and modeled extensively in the literature 
(see Reip urth et al. 19961 : Alencar et al.l 200ol : 
Folha et al.. .20011 and references therein). Ini- 
tial claims determined Balmer emission such 
as Ha to be attributed to outf lowing mate- 
rial (e.g., Hartmann et al. 1990l ) while, sub- 
sequently, magnetospheric accretions models 
have inve stigated their origin in terms of in- 



fall (e.g., lAlencar et alll2000l l. Indeed, the 
higher excitation Pa/? and Br7 lines have 
been found to correlate highly with accre- 
tion luminosity in classical T Tauri st ars and 



brown dwarfs (IMuzerolle et al.1 Il998al 1998b; 



Natta et al.l |200J). As a result, they have 



been used as mass accretion flux indicators in 
cases of highly embedded Class 0/1 sources 
where extinction is high, thus rendering the 
UV excess flux (i.e., the usual accretion indi- 
cator) difficult to measure accurately (e.g.. 



Garcia Lopez et al. 20061 ). However, there 
has never been satisfactory agreement be- 
tween models and observations, or full con- 
sensus within the community as to the true 
origin of these lines. More recent studies have 
again suggested an origin in outflowing mate- 
rial. There is some evidence th at Pa/3 seems 



to tr a ce T Tauri je t acti vity (iTakami et al 



200ll : IWhelan et al.l liooi : 2009). Also 



m- 



terferometric observations imply the Br7 line 
of Herbig Ae/Be stars may come from the 
i nner disk region and/or an ejected wind 
(|Kraus et al.l 120081 '). Meanwhile, Br7 detec- 
tions within 100 AU of T Tauri stars have 
been suggested as originating in the outflow 



or sca ttered from an outflow cavitv lBeck et al 
(I2OO9I ). 



In this paper, we examine the bipolar jet 
for T Tauri star ThA 15-28 (hereafter Th 
28) near the jet base (<1" along the jet, 
i.e., within 170 AU of the star). This is 
close to the region where the plasma is ac- 
celerated and collimated (i.e., within a few 
tens of AU of the star) just after jet launch- 
ing takes pl ace. The origin al paper on this 
star was bv iKrautteil (|l986l l. hence Th28 is 
now also known as Krautter's star. It is 
a classical T Tauri star located in the Lu- 
pus 3 cloud ( The 1962), which is at a dis- 
tance of 170 pc lEggen et al.l (119831'). and is sur - 
rounded in nebulosity ( Graham et al. 19881 ). 
The star's systemic heliocentric velocity, mea- 
sured from photospheric lines, w as found to be 
-hSkms"^ (jGraham et"ami988l ). The outflow 
inclination angle with respect to the plane of 
the sky was measu red as 10°, and the jet po- 
sition angle is 98° (|Krauttei] Il986l l. A study 
of the gas physics, based on ratios of opti- 
cal forbidden emission lines, at 2" along the 
flow showed the recedi ng jet to be of h i gh ion - 
isation (i.e., a;e~0.5) Bacciotti et al. ( 19991 ). 
Recent optical images confirm high proper 
motion of the bipolar flow HH 228 confirm- 
ing that it is excited by Th 28, see Figured] 
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Wang et alJ (|2009l ). 

The fact that this system is close to edge- 
on makes it ideal for examining the initial por- 
tion of the jet. Indeed, this bipolar jet was in- 
cluded in a survey on T Tauri jet targets which 
reported indications of jet rotation based on 
Hubble Space Telescope Imaging Spectrograph 
(HST/STIS) spectra at optical a nd ultraviolet 
wavelengths (ICoffev et alJ l2004l : 2007). The 
jet base was also the subject of a diagnos- 
tic study of the gas physics, which accessed 
jet plasma conditions within the first arcsec- 
ond from th e star at high spatial resolution 
Coffev et al. ( 20081 ). These studies had im- 
portant implications for the role of jets as an- 
gular momentum extractors, allowing accre- 
tion to proceed until the star reaches its final 
pre- main sequence mass. 

To build a multi- wavelength picture of this 
bipolar jet close to the source, we now choose 
to examine the Th 28 jet in near-IR emis- 
sion. The near-IR is expected to access the 
atomic jet via [Fe II] lines and the molec- 
ular flow via H2 emission. The aim of the 
observations was to understand the physics 
that can be extracted from the diagnostics of 
the infrared lines. For example, the electron 
density of the jet close to the star usually 
saturates the optical [S II]6731/6716 A ratio, 
whereas the near-IR allows access to higher 
electron densities before saturation of the 
[Fe Il]1.64/1.60 fj,m ratio, since [Fe II] has 
a higher critical density ('^6x 10^ cm~^) than 
[S II] (~2xl0'* cm~^). Higher densities de- 
termined in the initial jet channel have im- 
plications for mass and momentum fluxes, 
and hence the underlying mass loading at 
the jet base. Second, it is important to un- 
derstand the relationship between the atomic 
and molecular emission, both of which are 
accessed via near-IR observations. Finally, 
observations with the perpendicular slit were 
intended to determine if rotation signatures 
are confirmed for the atomic and molecular 



flow in the near-IR, as have been reported in 
optical lines, since it is important to establish 
that internal jet kinematics are consistently 
observed across all jet emission lines. 

2. Observations 

Spectroscopic observations were conducted 
using the Very Large Telescope (VLT) UTl 
located at Paranal, Chile under observing pro- 
gram 077.C-0818A. The ISAAC instrument 
was employed to cicc[iiir6 dcitci cit near-infrared 
wavelengths from the bipolar outflow from 
the T Tauri star, Th28. The VLT/ISAAC 
spectrograph in the MR configuration with a 
120x0.3 arcsec^ slit has a spatial sampling of 
0'.'147 and a spectral resolution of 8 900-10 500 
at the investigated wavelength region. The 
latter translates to a velocity resolution of 29- 
33kms~^. Where a line profile is intrinsically 
Gaussian, fitting allows an effective resolution 
which improves with climbing signal-to-noise, 
typically reaching one tenth of the resolution. 
Meanwhile, measurement of OH lines yield a 
deviation in absolute velocity calibration of 
~1.5kms~^ with respect to the theoretical 
value. 

Data were obtained in service mode over 
three nights: 2006 April 14th, 19th and 20th. 
A log of the observations is given in TableHl 
Observations were conducted with a single slit 
placed parallel to the bipolar jet using three 
filters, with a seeing of ~0'.'5 (FWHM) during 
the observations. The K filter was centered 
on 2.12 //m in order to capture H2 emission 
at 2.12 /um (i.e., at vacuum rest wavelength 
of 21 218.36 A). The i7 filter was centered on 
1.62 //m in order to capture [Fe II] emission 
at both 1.60 ^m (i.e., at vacuum rest wave- 
length of 15 999.086 A) and 1.64 //m (i.e., at 
vacuum rest wavelength of 16 439.98 A. The 
J filter was centered on 1.29 /^m in order to 
capture Pashen Beta emission (i.e., at vac- 
uum rest wavelength of 12 821.610 A). Obser- 
vations were also conducted with a single slit 
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placed perpendicular to each lobe of the bipo- 
lar jet at 0'.'5 (i.e., 86 AU deprojected dis- 
tance) from the star itself. An offset of Cf.'5 
was requested, based on the instrument PSF 
and assuming a seeing of (/.'G. This arrange- 
ment was chosen to provide a clear spectrum 
transverse to the jet axis and close to the star 
but without contamination from star light. 
However, the actual seeing was larger. The 
receding lobe was observed in the H and K 
filters, with a seeing on the night of C'T and 
C/.'8 respectively, while the approaching jet was 
observed in the H filter with seeing of C'T. 
We discuss the consequences of this in Sec- 
tion[3l The standard AB B'A' nodding tech- 
nique was employed during observations to al- 
low good sky subtraction while maximising 
available observing time. Using the optimum 
exposure time of 300 s, a total of eight individ- 
ual exposures were made in each filter, and co- 
added to achieve the required signal-to-noise. 

Data reduction was carried out using IRAF 
tools, and included cosmic ray rejection, flat- 
field correction using the lamp exposures, 
wavelength calibration using the telluric OH 
lines, distortion rectification using a star trace 
from the ESO archive, double sky subtrac- 
tion (or AB subtraction) using the nodded 
exposures, telluric correction and flux cal- 
ibration using the spectroscopic standard. 
A heliocentric radial velocity correction of 
-|-20kms~^ for all observations was deter- 
mined using IRAF tools. This was applied to 
all subsequent radial velocity measurements 
along with a correction of -5 ki ns~^ to negate 



the st ar's heliocentric velocity iGraham et al 



(|l988l ). 



Finally, we also present previously unpub- 
lished HS T/STIS Ha data along side the other 
H I lines. These observations were conducted 
on 2002 June 22 (proposal ID 9435), with 
the slit placed perpendicular to each jet lobe 
at O'.'S (52 AU deprojected distance) from the 
star. An aperture of 52x0.1 arcsec^ was used 



with the G750M grating which gave a radial 
velocity sampling of 25kms~^ and diffrac- 
tion limited spatial resolution of O'.'l. Spec- 
tra were obtained of the blue- and red-shifted 
lobes, using exposure times of 2200 and 2700 
s, respectively. They included several forbid- 
den lines (which w ere previously published; 
Coffev et al.l (|2004l )) along with Ha which we 
present here. Data were calibrated through 
the standard HST pipeline, subtraction of 
the stellar continuum was performed, and 
hot /dark pixels were removed. 

3. Results 

The Herbig-Haro outflow from Th28, 
named HH 228, consists of a bipolar microjet 
close to the star, two bowshocks 39" and 45" 
either side of the star aligned close to east- 
west, and a third fainter bowshock farther to 
the east at 95" from the star. The overall posi- 
tion angle of the flow is almost east-west, with 
the brighter red lobe of the microjet traveling 
towards the west, as illustrated in Figured) 
The jet lies close to the plane of the sky, and 
consequently has low radial velocities of a few 
tens of kms~^. High proper motions of 250 
- 350 kms^^ co nfirm that these o bjects are 
driven by Th 28 IWang et all (|2009l l. 

For the first time, we have detected the 
bipolar microjet in near-IR emission via spec- 
tra taken with the slit both parallel and per- 
pendicular to the flow direction. We will first 
describe the 'parallel' configuration spectra, 
which illustrate the global morphological and 
kinematic properties of the flow close to the 
source, and then we will focus on the 'perpen- 
dicular' spectra, which highlight further im- 
portant peculiarities. 

3.1. Parallel-slit Configuration 

In the parallel case, the position-velocity 
(PV) plots, collected in Figure[21 show detec- 
tions in [Fe II], Br7, Pa/3 and H2. In these 
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Fig. 1. — Continuum subtracted [S II] image of HH 228, adapted from IWang et alJ (|2009l ). with 
our VLT/ISAAC slit positions overlaid. 
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Slit position Filter Central wavelength Spectral Domain Resolving Power Seeing Exposure time 

(/im) (m™) (arcsec) (s) 

Parallel to bipolar jet K 232 2.06 - 2.18 8 900 05 2400 

H 1.62 1.58 - 1.66 10 000 0.6 2400 

J 1.29 1.26 - 1.32 10 500 0.5 2400 

Perpendicular to red lobe K 2.12 2.06 - 2.18 8 900 0.8 2400 

H 1.62 1.58 - 1.66 10 000 0.7 2400 

Perpendicular to blue lobe H L62 1.58 - 1.66 10 000 07 2400 



Table 1: Details of observations of the jet from T Tauri star ThA 15-28 (a.k.a. Th 28) conducted 
with VLT/ISAAC in service mode over 3 nights in April 2006. 




-300 -20O -IW l(ffl 200 300 -300 -200 -100 100 200 300 -300 -200 -100 100 20O 300 -300 -200 -100 100 2(» 300 -300 -200 -100 100 200 300 

Radial Velocity (km s"') Radial Velocity (km s'') Radial Velocity (km s'') Radial Velocity (km s'') Radial Vekicity (km s'') 



Fig. 2. — PV plots of the Th 28 bipolar jet with the long slit parallel to the jet axis, illustrating 
the extent of emission with distance from the star. Each plot has ten linearly-spaced contour, with 
contour levels from 3a to 70, 890, 240, 70, 670c7 respectively. 
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plots, the position of the central star was de- 
termined from a Gaussian fit to the contin- 
uum, which was then subtracted. 

The various emission lines all show a wide 
emission peak around the star, but since they 
probe different physical conditions, they ex- 
tend to different distances from the source. 
In the atomic lines, the red lobe is brighter, 
and corresponds to negative distances along 
the jet in Figure[2j [Fe II] traces both lobes of 
the bipolar jet, with [Fe II] 1.64 ^m emission 
from the red lobe being stronger and persist- 
ing along the jet to 6", compared to 2" for the 
blue lobe. The Pa/3 emission clearly extends 
to 3" while Br7, which is intrinsically fainter 
than Pa/3, is closer to the jet base and we do 
not see it beyond 1.5" along the receding flow. 
The fact that Pa/3 and Br7 are clearly seen 
at a distance from the source tracing the red 
lobe in the PV diagrams is very interesting, 
as these lines have been often commonly at- 
tributed to accretion phenomena. Meanwhile, 
the H2 emission presents a different behaviour 
with respect to the atomic lines: it is more 
confined to a narrow range of low velocities, 
and is more extended on the side of the blue 
lobe. 

To further examine the structure and kine- 
matics indicated by these lines, the flux was 
spatially binned to match the seeing. The spa- 
tial profile wings extend beyond a simple spa- 
tial Gaussian fit of the line emission (dashed 
curve) implying origin in an extended source, 
FigureO left column. Corresponding velocity 
profiles along the jet were determined from 
Gaussian fits in the velocity direction per- 
formed at different distances from the source, 
FigurelU right column. We see that [Fe II] 
travels with a velocity of about -|-30kms"^ in 
the red lobe and -50kms~^ in the blue lobe. 
The plots also show that Br7 and Pa/3 travel 
with similar velocity to the well known atomic 
jet tracer [Fe II]. These permitted H lines ap- 
pear to be produced in the jet to distances of 



zbO'.'G from the star for Br7, and from -3" to 
+1'!2 from the star for Pa/3. In the same CCD 
frame with Pa/3, four other iron lines of simi- 
lar excitation to [Fe II] 1.64 //m, but lesser in- 
tensity, were also observed in the parallel slit, 
namely [Fe II] 12706.94, 12791.255, 12946.232 
and 12981.279 A( vacuum wavelengths). All 
of these show velocity profiles consistent with 
that of the [Fe II] 1.64 ^m line. More signif- 
icantly, the higher excitation [Fe II] 2.13/Lim 
line was also detected. This is an important 
confirmation that the similarly high excitation 
hydrogen lines, Br7 and Pa/3, originate in the 
jet itself. 

Finally, we find that the H2 emission unex- 
pectedly shows blue-shifted velocities of about 
-lOkms"^ in the direction of both the reced- 
ing and approaching jet lobes. Its spatial dis- 
tribution is also different to the atomic lines, 
in that it is more evenly distributed about the 
system, and slightly more extended toward 
the blue lobe. In addition, the H2 line is nar- 
rower than the atomic lines, indicating that 
the emitting material here is colder and/or 
less affected by turbulent motions. 

3.2. Perpendicular-slit Configuration 

To understand in further detail the condi- 
tions close to the jet base, we examine each 
jet lobe via our data taken with the slit per- 
pendicular to the jet axis at a distance of (f.'5 
from the source. In Figure|l]we first show flux- 
calibrated one-dimensional extracted spectra, 
that globally present the detections in near- 
IR molecular and atomic lines. We find H2, 
[Fe II] lines and again, interestingly, Br^. 
While both lobes were detected in [Fe II] lines, 
emission was fainter for the approaching jet 
than for the receding jet, highlighting the fact 
that the former is more embedded. Note that 
the spectral FWHM of the H2 line in the re- 
ceding lobe is much smaller than that of the 
other lines, in agreement with the findings for 
the 'parallel' spectra. Perpendicular slit H2 



7 



Th 28jet-[Fell]1.64nm 



fwhm continuum 0.6" 
fwhm line 0.9" 




Distance from star (arcsec) 
Th 28 jet - Pap 1 .28|im 



Iwlim continuum 0.5" 
fwhm line 0.6" 



-5 5 

Distance from star (arcsec) 

Th 28 jet-BrY2.16|im 



fwhm continuum 0.5" 
fwhm line 0.6" 



fwhm continuum 0.5" 






fwhm line 0.6" 








L 


,_ — ^ » -i 



-5 5 

Distance from star (arcsec) 

Th28jet- H 2.12M.m 



Distance from star (arcsec) 



Th28jet-[Fe II] 1.64nm 


* 








1 

_L 

1 
1 




W {red lobe) 


1 

1 $ 

1 , , , , 


E (blue lobe) 


-5 




Distance from star (arcsec) 
Th 28 jet - Pap 1 .28|im 


5 




1 
1 

\J 

1 




W (red lobe) 


1 
1 
1 
1 

1 , , , , 


E (blue lobe) 


-5 




Distance from star (arcsec) 
Th28jet-Bi72.16nm 


5 




o 




W (red lobe) 




E (blue lobe) 


-5 




Distance from star (arcsec) 
Th 28jet- H^2.12M.m 


5 




1 
1 
1 
1 
1 
1 
1 

^-^^^ 




W (red lobe) 


1 , , , 


E (blue lobe) 


-5 




Distance from star (arcsec) 


6 



Fig. 3. — Spatial intensity and velocity profiles of Th 28 bipolar jet in near-IR lines through a 
slit parallel to the flow axis. The wings of the profile extend outside a simple Gaussian fit to the 
emission (dashed line) indicating an extended source. Flux is spatially binned to match the seeing. 
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observations were not conducted for the ap- 
proaching jet. 

The main emission properties derivable 
form extracted spectra are collected in Ta- 
ble[21 where for each line we list the calibrated 
flux and signal-to-noise, the velocity of the 
emission peak and the velocity FWHM, and 
finally the jet spatial width. The jet width 
was resolved only in H2 (FWHM of 1'.'6) under 
the seeing during the observations of not bet- 
ter than ff.'7. The fact that the atomic near- 
IR lines were spatially unresolved precluded 
any conclusions on jet rotation for these lines, 
since this would be determined from radial ve- 
locity gradients detected across the jet width. 
Nevertheless, the available near-IR perpen- 
dicular spectra were very useful in identifying 
relevant physical conditions in the jet. 

The information of Table[2] is further com- 
plemented by FigureO where we present the 
PV diagrams of [Fe II] 1.64 ^m line in both 
lobes and Br^ in the red lobe. For com- 
parison with optical emission, in Figure[6] we 
also present a PV plot of the Ha emission 
in HST/STIS spectra, taken with a similar 
instrument configuration, but as yet unpub- 
lished. 

The PV diagrams reveal a clear velocity 
asymmetry between the red and blue lobe. 
The high-velocity wing of the iron line in the 
blue lobe reaches -200kms~^ while the red 
lobe reaches +120kms~^. The asymmetry is 
reflected in Table[2] which reports the peak ve- 
locities via Gaussian fits to the integrated ve- 
locity profiles, and yields +25kms^^ for the 
receding jet and -15kms~^ for the approach- 
ing jet [Fe II] emission. A similar behaviour is 
observed in the optical lines measured with 
HST/STIS, FigureEl In this case too, the 
emission on the blue side is fainter but it ex- 
tends to higher velocities, with limiting values 
similar to the iron lines. The peak velocities 
and resulting asymmetry for Ha are greater 
than the near-IR lines however, with values 



measured as -127 km s ^ and +62kms ^ for 
the blue and red lobes respectively. 

On closer examination of the PV plots, it 
is clear that near-IR atomic emission covers a 
broad velocity range either side of zero. This 
is confirmed by a comparison of the peak ve- 
locity with the FWHM values of Tabled a 
spread which is too large to be explained by 
the velocity resolution. Indeed, the Ha PV 
plot shows that in the blue lobe the emis- 
sion appears to be double-peaked, with two 
distinct velocity components, one blueshifted 
and one redshifted. A double Gaussian fit to 
the profile reveals velocities consistent with 
the red and blue lobe velocities. It seems that 
the slit positioned over the blue lobe also takes 
in reflected light from the red lobe. This is 
further confirmed by the spatial asymmetry 
of the redshifted peak which matches exactly 
with the obvious spatial asymmetry of the re d 
lobe PV plots in [O I] ICoffev et al.l too4 ). 
and is also coincident with an electron den- 



sity a symmetry across the axis ICoff^ev et al.l 
(|2008l ). A similar transverse spatial asym- 
metry could not be identified in the bright 
[Fe II] 1.64 /im line, due to overlap with a 
night sky line. The scattered light is likely 
to be caused by the nebulosity around the 
blue lobe. Note that there is a less obvious 
refiection in the slit of the red lobe, which 
is compatible with the blue lobe being more 
embedded. This refiection from the opposing 
lobe was not encountered in the other targets 
examined under these HST proposals with a 
perpendicular slit configuration (proposal ID 
9435; 9807). The implications for the parallel 
slit spectra are simply that any contamina- 
tion from the opposing lobe will result in an 
under-estimation of the radial velocity values 
and an overestimation of the fiux, but there 
are no implications for the interpretation. In- 
deed, the main concern here is the true ori- 
gin of Br^ in the perpendicular-slit "jet" spec- 
trum. Our slit was positioned at C/.'S from the 
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Fig. 4. — Extracted one-dimensional spectra of near-IR atomic and molecular lines, derived from 
jet spectra with the slit perpendicular to jet propogation. The spectra reveal single-peaked emission 
in H2, [Fe II] and Br^ originating in the region of the microjet. The H2 emission is stronger but 
much narrower, implying origin in a medium colder than the atomic lines. 



Target Emission Flux FWHM S/N v Seeing Jet width 

(^m) (10^"'"® ergcm"^ s"-"-) (kms"-*-) (kms^-*-) (arcsec) (arcsec) 

Red lobe H2 A2.12 7.1 ±0.1 45 112 OS L6 

[Fe II] A2.13 1.0 ±0.3 138 6 +23 0.8 unresolved 

Br7A2.16 3.1 ±0.3 187 11 +6 0.8 

[Fe II] A1.64 82.5 ±0.3 134 315 +24 0.7 

[Fe II] A1.60 16.7 ±0.3 139 69 +18 0.7 

Blue lobe [Fe II] A1.64 17.1 ±0.3 215 49 -16 0.7 

[Fe II] A1.60 3.1 ±0.3 215 12 -13 OT " 



Table 2: Intensity and kinematical information derived from near-IR emission from the Th 28 bipolar 
jet obtained with a slit perpendicular to the flow direction at C/.'S from the star. The jet width is 
spatially unresolved in all atomic lines, contrary to the H2 emission. Fluxes may be over-estimates, 
and velocities under-estimates due to scattered light (see Section [3^ . 
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Fig. 5. — PV plots of the Th 28 bipolar jet with the long slit perpendicular to the jet axis. Each 
plot has ten linearly-spaced contour, with contour levels from 3 to 60, 180 and 2Qa respectively. A 
night sky line disrupts the [Fe II] emission. 




Fig. 6. — Th 28 bipolar jet observed in Ha via HST/STIS at O'.'S from the star, with a long slit 
positioned perpendicular to the flow direction. Contours in units of ergcm"^ s~^ arcsec^^ are 
from 3x10-15 to 1.22x10-^^ and 6.0x10-^^ respectively. 
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star in each lobe, to allow for the instrument 
PSF and to avoid stellar flux. Unfortunately, 
since the seeing on the night was in fact C/.'? 
- (X'S, the spectra unavoidably include contri- 
butions aside from the jet itself. As expected, 
stellar continuum emission is seen in these 
spectra, but was subtracted prior to analy- 
sis of the lines. Meanwhile, we overcome the 
difficulty of scattered light by relying on ve- 
locity measurements to indicate whether the 
emission traces the jet. The Br^ velocity of 
-f6kms~^ in the receding jet, with a velocity 
error of 1.5 kms~^, combined with a detection 
of [Fe II] 2.13/im emission in the same region, 
strongly argues that the Br^ line truly trace 
outflowing gas. 

Finally, we consider molecular hydrogen 
detected in the perpendicular slit spectrum. 
Figure[7] shows a PV diagram, along with 
the spatial profile and the velocity profile ob- 
tained by Gaussian fitting. As we saw in 
the parallel-slit spectra, the characteristics of 
the H2 emission are unlike the ionic lines. 
It is broad enough to be spatially resolved, 
it exhibits a narrow spectral profile extend- 
ing over 45kms~^, and it supports the unex- 
pected finding of the parallel-slit data that the 
emission traces a blue-shifted gas (-12 kms~^) 
in the region of the receding jet. From these 
data, the H2 gas detected around the base of 
the bipolar jet is entirely blueshifted to about 
-lOkms"^ in both the parallel and perpen- 
dicular slit data regardless of whether the slit 
was positioned over the approaching or reced- 
ing jet. It appears that only the blue lobe is 
emitting in H2 , and that we are detecting scat- 
tered light from the blue lobe when observing 
in the direction of the red lobe. This points 
to an asymmetry in extinction and/or excita- 
tion between the two lobes. Since the perpen- 
dicular slit emission is spatially resolved, we 
may extract a velocity profile across the gas 
to determine if there are possible signs of ro- 
tation. Figure[71 right panel, shows a slight 



systematic gradient in radial velocity, which 
might give a borderline indication of rotation 
in the H2 gas. However, the sense of implied 
rotation is opposite to that measured for op- 
tical atomic jet lines resolved with HST/STIS 
which presented skews in the profile of 10- 
25km s~^ for the same target Coff^ev et al.l 
(I2OO4I I. However, given the nature of this 
emission, any kinematic asymmetry is more 
likely caused by environmental factors, rather 
than intrinsically related to the jet launching 
mechanism. 

4. Discussion 

4.1. Multi- wavelength jet picture 

These near-IR spectra complement HST/STIS 
optical and near-UV observations for the same 
target, conducted with the same slit position 
and orientation (proposal ID 9435; 9807). All 
three datasets agree that atomic emission is 
fainter for t he approaching jet than for the 
receding jet (jCoffev et al.1 l2004l : 2007), high- 
lighting the fact that the former is more em- 
bedded. 

Examining the combined topology, we see 
that the width of the receding jet was de- 
termined to be O'.'S and 0^.^ 4 for optical [O Ij 
and [S II] lines respectivelv ICofFev et al.ll2004l . 
while in the near-UV M g II lines it was 
0'.'l-0'.'2 for both jet lobes Icoffev et al.l 12007 . 
Meanwhile, the near-IR atomic lines are not 
spatially resolved with a seeing of 0'.'7-0'.'8, 
and the molecular gas presents an outflow 
width of 1'.'6. This multi-wavelength pic- 
ture illustrates that the jet constitutes sev- 
eral layers of nested material with the near- 
UV emitting region at the centre, the opti- 
cal and near-IR emitting material further out 
and flnally a molecular layer. Indeed, this 
structure is commo n among jets ( e.g., t he jet 
from DG Ta u . see ICoffev et al.l ^200l\ ) and 
Takami et all (j20oi))- 

Velocities are consistent across all wave- 
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TH 28 red lobe . 
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Fig. 7. — Left: PV plot of H2 in the direction of the Th 28 receding jet (via ten hnear contours from 
3 to 104(T. ). Middle: the spatial profile of the molecular emission shows a slight asymmetry in 
emission with respect to the flow axis. Right: transverse velocity profile (spatially binned according 
to the seeing) of the resolved molecular emission, indicating clearly that the emission is entirely 
blueshifted even though the slit is positioned over the receding jet lobe. 



lengths regimes. We have seen that [Fe II] 
travels at about +30kms~^ in the red lobe 
and -50 km s~^ in the blue lobe. This asymme- 
try is in excellent agreement with the asym- 
metry of +30 and - 60 km s~^ reported for opti- 
cal emission lines dCoffev et al.ll2004l ^. Mean- 
while, the Ha emission is traveling with peak 
velocities of -127kms~^ and +62kms~^ for 
the blue and red lobes respectively. (Recall 
that velocity measurements may be under- 
estimated, depending on the contribution of 
scattered light in each emission line.) Such 
velocity asymmetries between the lobes of a 
bipolar jet are common, and have been re- 
ported in many other ca. ses, e.g., the RW 
Au r ietlWoitas et al.1 (I2OO2I I and the LkHa233 
jet 



Melnikov et al 



(|2008l ). These asymme- 
tries are still not fully explained, but recent 
estimates of the corresponding mass outflow 
rates, which transpire to be similar in the two 
lobes, appear to indicate that the asymme- 
try is due to environmental conditions on the 
two sides of the system rather than to an in- 
trinsic asymme t ry of the accelerating engine 
Melnikov etHI (|2009l '). 



base of the bipolar jet is entirely blueshifted 
to about -lOkms"^ in both the parallel and 
perpendicular slit data regardless of whether 
the slit was positioned over the approaching 
or receding jet. It appears that only the 
blue lobe is emitting in H2, and that we are 
detecting scattered light from the blue lobe 
when observing in the direction of the red 
lobe. This points to an asymmetry in extinc- 
tion and/or excitation between the two lobes. 
Also, the receding jet is more luminous in 
atomic lines. This situation may arise because 
the approaching jet is more embedded, or be- 
cause the receding jet is somehow intrinsically 
more ionised. However the near-UV Mg II 
doublet line ratio is one in the blue lobe (com- 
pared to two in the red lobe), which implies a 
higher column density of Mg II. This optical 
thickness in the blue lobe indicates that it is 



The H2 gas which we detect around the 



more embedded ICoffey et al.l (120071 ). and thus 
supports an asymmetric extinction. 

4.2. Electron density and mass flux 

The ratio of near infrared iron lines allows 
us to derive the density in the region, which 
is usually high enough to saturate the com- 



13 



monly used [S II] doublet at optical wave- 
lengths. Figure[8] gives the electron density, 
Ue, derived from [Fe Il]1.64/1.60^m ratio, as- 
suming a temperature 10^ K. The upper panel 
gives the normalised intensity of the parallel 
slit observations, and the lower panel gives 
electron densities derived from parallel slit 
fluxes which were binned to match the see- 
ing. Over-plotted are electron densities de- 
rived from the perpendicular slit data for the 
blue and red lobe. Measurements along the 
jet show that close to the star the electron 
density saturates the [Fe Il]1.64/1.60//m ratio, 
indicating a lower limit of 610^cm~^. How- 
ever, this value rapidly falls off within the 
first arcsecond. At a distance from the star 
of C/.'S, the electron density is already down to 
2 (±0.1) 10^ cm~^ in the red lobe, as measured 
in both parallel and perpendicular slit data. 
Within errors, both lobes can be said to have 
the same density, although a density asymme- 
try may exist given the errorbars for the blue 
lobe i.e., ne=1.4(±0.4) 10*^ cm"^. Meanwhile, 
the electron density profile across the jet is 
unresolved, and so could not be compared to 
the density asymmetry derived from the [S II] 
doublet in the optical rang e from HST/STIS 
spectra ICoffev et al.l ( 20081 ). The latter report 
a typical electron density for the red lobe at 
O'.'S from the star of 1 10^ cm~^, although this 
does not take into account that one side of the 
jet shows a region which saturates the [S II] 
doublet (i.e., lower limit on electron density of 
210^cm~'^). Finally, the higher electron den- 
sity closer to the star derived from the iron 
lines is consistent with the fact that we ob- 
serve the higher excitation H I lines of Pa/3 
and Br7 in this region, as high densities are 
required to populate the upper levels of these 
permitted transitions. 

Mass flux estimates are of significant im- 
portance in order to access, for example, 
the drive of molecular flows, their contribu- 
tion to angular momentum extraction, and 




-2-10 1 
Distance along jet (arcsec) 

Fig. 8. — Top: Spatial profile of the nor- 
malised flux from [Fe II] 1.64 and 1.60/im 
emission (solid and dashed lines respectively) 
from parallel slit observations. Bottom: Elec- 
tron density derived from the [Fe II] 1.64/1.60 
line ratio as a function distance from the 
source, using parallel and perpendicular slit 
datasets (diamonds and crosses respectively). 
For the parallel slit, the flux was spatially 
binned to match the seeing. Here, the elec- 
tron density could not be derived for the blue 
lobe (positive distances along the jet) due to 
insufficient intensity, unlike the perpendicular 
slit data which takes in flux from the full jet 
width. 
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the impact of jets and outflows on their 
environment. Furthermore, accurate esti- 
mates of the mass flux rate close to the ejec- 
tion point of the jet have important impli- 
cations for the underlying mass loading at 
the jet base. We determine the mass flux 



rate, Mjet, following method B of lNisini et al 



(|2005l l and IPodio etaTI (|2006l ). The method 



derives Mjet from the observed luminos- 
ity, L{line), of optically thin lines which 
is proportional to the mass of the emit- 
ting gas: Mjct = fJ-mn {nHV)vt/lt, with 

uhV = L{line) (jiuAifi^^^ , where 
/i=1.41 is the mean atomic weight, mn is the 
proton mass, uh the hydrogen density, V is 
the volume filled by the emitting gas, vt and 
It are the tangential velocity and length of the 
knot, Ai and fi are the radiative rate and the 
upper level population relative to the consid- 
ered transition, and ^ and ^ are the ion- 
isation fraction and the relative abundance 
of the considered species. This method is 
affected by uncertainties in absolute calibra- 
tions, extinction, and distance. On the other 
hand, it does implicitly take into account the 
volume filling factor. Alternatively, Mjet may 
be derived from estimates of the total density 
and jet radius. Here, the jet section is as- 
sumed unifor mly filled at a give n density (see 



Method A of IPodio et all ^2Q0(i )). which may 



result in an upper limit to the real mass flux. 

To estimate Mjet at the base of the jet from 
Th 28, we consider the [Fe II] 1.64 //m line flux 
observed using the perpendicular slit. Since 
the red lobe is the brightest, we can measure 
line flux with high signal-to-noise and thus ob- 
tain an accurate estimate of electron density. 
We have no estimate of the visual extinction 
toward Th 28, so we assume Ay=0, i.e., we 
do not correct the derived luminosity for pos- 
sible reddening effects. Finally, we assume 
a temperature, T, of 10^ K to compute the 
upper level population (fi), we assume iron 
to be completely ionised (^=1) and we as- 



sume a solar abundance of iron 



Asplund et al 



2.82 10"^ 
Based on these as- 
we find a value of ~2 10~^ and 1 



(|200i)) 

sumption 

10~^Mo yr~^ for the mass flux rate of the re- 
ceding and approaching jet respectively. The 
red lobe yields twice the mass flux of the blue 
lobe. This is due to the lower apparent flux 
of the blue lobe (~20% of the red lobe flux), 
which is only slightly compensated by the fact 
that the velocity is higher in the blue lobe (i.e., 
twice the red lobe velocity). The mass flux 
asymmetry however does not account for an 
asymmetric extinction between the two lobes. 

Our values are lower than the typical mass 
flux rate o f ^10~^ MfDyr"^ rep orted for T 
Tauri stars iHartigan et al.l (jl995l ). The out- 



come is a lower limit on the mass flux rate for 
two reasons. First, we do not have an esti- 
mate for the extinction and, second, the iron 
gas phase abundance may be d epleted with 



respec t to the solar abundance (jPodio et al 



(I2OO6I ): 200 8). The valu e repo rted from op- 



Coffev et al.l (|2008l l is 1.2 10" 



tical lines 

M0yr~^, calculated by using the total den- 
sity and jet radius, i.e., method A. This order 
of magnitude difference in mass flux rate be- 
tween the two methods (i.e., relying on the 
luminosity as opposed to electron density and 
ionisation determinations) is the typical dif- 
ference found in previous studies which com- 
pare the two method. Note that for this jet 
target, the flux of the optical lines is roughly 
two orders of magnitude higher than the 
[Fe II] 1.64;um flux. This can be due to iron 
depletion and/or intrinsically fainter emission 
in the iron lines with respect to the optical 
lines in this target. 

4.3. H I Origin 

There is ongoing debate as to whether the 
origin of hydrogen line emission close to the 
star is accretion and/or ejection related (Sec- 
tion[T]). Our results demonstrate that a com- 
ponent of both Pa/? and Br7 emission is un- 
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doubtedly coming from the outflow. This is 
evident from the clear spatial extension seen 
in the parallel slit data, and from the veloc- 
ities found in these spatially extended wings. 
The finding is interesting given that these are 
high-excitation lines. Both lines present simi- 
lar velocities to the [Fe II] lines which we know 
to trace the jet, and are further confirmed to 
originate in jet by the detection of the sim- 
ilarly high excitation [Fe II] 2.13 //m line. 
This is a highly excited iron line, which has 
been detected in jet knots in conjunction with 
Br7 ( see e.g., HH 34 in iGarcia Lopez et al 



(2003)). It is a further indication that Br7 
originates in the jet. The velocities and fluxes 
of these H I lines measured close to the star 
(Table[2]) provide important constraints on 
models. The proportion of flux originating in 
the outflow could not be determined because 
these lines arise in the higher density jet re- 
gion close to the star which suffers terms of 
seeing limitations and scattered light. Finally, 
the importance of identifying Pa/3 and Br7 as 
originating in the outflow lies in the fact that 
they are permitted transitions in the near-IR. 
Permitted transitions emit in high density re- 
gions where forbidden lines are quenched, and 
near-IR lines are less affected by extinction 
than optical lines. Hence, the near-IR lines 
Pa/3 and Br7 become very important probes 
of the embedded central engine of protostars. 
The development of high resolution interfer- 
ometric instruments in the infrared, such as 
VLTI/ AMBER, means it is important to iden- 
tify in advance useful probes of the central 
engine. 

5. Conclusions 

Our study provides the flrst information on 
the near-IR spectrum of the microjet from Th 
28, revealing emission in [Fe II], Pa/3, Br7 
and H2. This near-IR dataset complements 
the optical and near-UV data taken for the 
same target with HST/STIS. The Th28 out- 



flow now constitutes the only bipolar jet sys- 
tem to date for which we will have a complete 
spectroscopic dataset sampling all the coaxial 
layers of the jet, and for both jet lobes. 

The molecular H2 emission was detected 
surrounding the base of both jet lobes, but 
its blueshifted velocity of -10 km s~^ indicates 
that only the blue shifted lobe is emitting 
while light is scattered in the direction of the 
red lobe, highlighting an asymmetric extinc- 
tion and/or excitation between the two lobes. 
Outflow radial velocities in ionic lines are in 
agreement across all wavelengths though the 
velocities, and the velocity asymmetry be- 
tween the lobes. The HST/STIS Ha emission, 
observed at 0'.'3 from the star, presents ra- 
dial velocities of -127kms~^ and +62kms~^ 
for the blue and red lobe respectively. Un- 
like optical and near-UV results, there is no 
indication of jet rotation in [Fe II] lines be- 
cause of insufficient angular resolution to re- 
solve the jet width. The iron ratio was found 
to saturate close to the star implying a lower 
limit on the electron densities of 6 x 10^ cm^ . 
However, the high density falls off rapidly to 
2 10^ cm^ within the first arcsecond along each 
flow, and leads to a value of 2 and 1 10~^ 
M0yr~^ for the mass flux rate of the re- 
ceding and approaching jet calculated from 
iron line apparent luminosity. This is an 
order of magnitude lower than that derived 
using optical lines, likely due to differences 
inherent in the method and possible deple- 
tion of iron in the gas phase. Finally, a ve- 
locity analysis, along with the detection of 
the high excitation [Fe II] 2.13 /xm line, re- 
veals that the jet is traced in Pa/3 and Br7 
within the flrst 150 AU from the star. These 
higher excitation hydrogen lines are usually 
attributed to accretion activity, but our ob- 
servations clearly demonstrate there is a com- 
ponent which traces outflow. These observa- 
tions, though seeing limited, have good spec- 
tral resolution and thus will complement the 
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high angular resolution of VLTI/ AMBER in- 

terferometric observations when sensitivities 
are sufficient to observe T Tauri jet emission. 
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